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ABSTRACT 

We present radio and X-ray observations of the nearby Type lib Supernova 2013df in NGC4414 
from 10 to 250 days after the explosion. The radio emission showed a peculiar soft-to-hard spectral 
evolution. We present a model in which inverse Compton cooling of synchrotron emitting electrons 
can account for the observed spectral and light curve evolution. A significant mass-loss rate, M ss 
8 x 10 -5 Mq /yr for a wind velocity of 10 krn/s, is estimated from the detailed modeling of radio and 
X-ray emission, which are primarily due to synchrotron and bremsstrahlung, respectively. We show 
that SN 2013df is similar to SN 1993J in various ways. The shock wave speed of SN 2013df was found 
to be average among the radio supernovae; v s h/c ~ 0.07. We did not find any significant deviation 
from smooth decline in the light curve of SN 2013df. One of the main results of our self-consistent 
multiband modeling is the significant deviation from energy equipartition between magnetic fields and 
relativistic electrons behind the shock. We estimate e e = 200es- In general for Type lib SNe, we find 
that the presence of bright optical cooling envelope emission is linked with free-free radio absorption 
and bright thermal X-ray emission. This finding suggests that more extended progenitors, similar to 
that of SN 2013df, suffer from substantial mass loss in the years before the supernova. 

Subject headings: radiation mechanisms: nonthermal - radio continuum: general - supernovae: general 
- supernovae: individual (SN 2013df, SN 1993J) 


1. INTRODUCTION 


With the recent advent of high cadence (< 1 day) op¬ 
tical transient searches, our observational understanding 
of the nature of stellar explosions has bloomed. Through 
these efforts, exotic breeds of supernovae have been re¬ 
vealed with peculiar light-curve and spectral properties, 


explosion (e.g., 

Kasliwal et al. 

2010 

Drout et al. 

2013). 

The rare and intriguing class o 

Type lib supernova (SN 


lib) showing evidence for both hydroge n and helium 


in ea r ly spectroscopic observa tions (Ensman & Woosley 
1987[ Filippenko 1988 1997) - has become one focus 


point of these new optical transient surveys. The lu¬ 
minosity and evolution of the early emission places strict 
constraints on the physical parameters (e.g., mass, ra- 


system (Baron et al. 

1993 Filippenko, Matheson, & Ho 

1993 

Swartz et al. 

L993) as it maps to the cooling en- 


breakout from the stellar surfa ce (|Katz, Budnik, & Wax- 
man||20i0 Nakar fe Sari||20l0 ). 


Here we present radio and X-ray observations of the 
recent type lib SN 2013df in NGC4414 to investigate 
the nature of its progenitor system and the pre-explosion 
mass loss history. Based on detailed modeling of the en- 
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tire dataset, we clearly show that the SN is essentially 
an analogue to SN 1993J and stems from an extended 
progenit or star, in line with th e results of pre-explosion 
imaging) Van Dyk et al. 2014|) and post-explosion UV 
observations (Ben-Ami et al. [2014). As in the case of 
SN 1993J, our early radio and X-ray observations point 
to dense material in the local environs as evidenced by 
inverse Compton processes. We show that these obser¬ 
vations imply diversity among the progenitors and may 
provide crucial diagnostics to distinguish between possi¬ 
ble sub-classes of Type lib SNe. 

SN lib were originally identified as a unique class 
of core-collapse supernovae two decades ago with the 
discovery and long-term m ulti -wave length stu d y of 
SN 1993J in M81 ( FUiiJijeiiko^^Mathesom&^HoJlOQd van 


Dyk et al.|1994 |Fooley &; Green|1993 " Filippenko|1997 ). 
There is general consensus that the progenitor star of 


general consensus tliat tlie progenitor star of 
SN 1993J was a massive K type super-giant star with 
R ~ 600 R ( 7) and M « 13—22 based on pre-explosion 


imaging (Aldering, Humphreys, & Richmond 

19941 

o 

o 

hen, Darling, & Porter 19<J3|; Maund &; Smartt 

2009; 

Fox 


ied in depth and the early optical emission was observed 
to decline over the first ~ 6 days prior to rebrightening 
due to radioactive decay products. In addition, the un¬ 
usual nature of the late-time spectral evolution proved 
indicative of circumstellar interaction with dense mate¬ 
rial ejected prior to stell ar collapse the origin of which 



normal H-rich Type II supernovae (SN II) and H-poor 
Type Ibc supernovae (SN Ibc) based on their spectro¬ 
scopic properties. 

One unifying feature and hallmark of the class of SN 
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lib is circumstellar interaction as evidenced by strong 
non-thermal emission in the radio and X-ray bands due 
to dynamical intera ction of th e shockwa ve w ith the_sur- 
rounding medium (jChevalier 1982 Chevalier & Frans- 
’20031). In addition to SIN 11J9.1 J. such signatures 


son 


have been observed for Type lib SN 200 1ig (Ryder, 


Murrowood, fc Stathakis ] 2006), SN 2003bg (Soder befgT 

a.il||2006|) , SN2008ax (IRoming 

ioderber 


et al. 

2009) 

et al. 

2012; 

jevic et al. 


iIN 2Ulldh (ISoaerberg et ai.| 
e_Witt et al.||2()15|). SN 2pl 


21)12[ :Kraiiss 

lei ([Milisavl- 


1>013), SN 2U12aul Gamble et al. 2014). In 


ot these cases, episodic liux density modula¬ 
tions of the non-thermal emission are observed and point 
to circumstellar medium (CSM) density fluctuations on 


a radial scale 5, 10 17 cm (Ryder et al. 


Soderberg, &; Chevalier| 2012[), possibly due to a binary 


2004 Wellons, 


companion or mass loss variability (due to stellar pulsa¬ 
tions) prior to explosion. 

Combining early optical observations with radio and 
X-ray data for a small sample of SN lib, |Chevalier fe| 


Soderberg (2010) proposed that the class of SN lib is di¬ 
verse and encompasses both “compact” (Type cllb; R ~ 
Rq) and “extended” (Type ellb; R ~ 100 — 1000 Rq) 
progenitor stars. The authors suggest that compact pro¬ 
genitors are distinct and can be identified based on the 
shockwave velocity as well as the luminosity and dura¬ 
tion of the optical cooling envelope phase. In this frame¬ 
work, SN 1993J was clearly identified as a SN ellb while 
most of the other Type lib’s were identified as being 
more closely aligned with a SN cllb classification, thereby 
bearing more resemblance to SN Ibc. We note that the 
binary classification into either SN cllb or SN ellb is 



2. OBSERVATIONS 

2.1. Radio Observations using Very Large Array 

SN 2013df was optically discovered on 7.87 June 2013 
by the Italian Supernovae Search Project, with an offset 
of about 32 , . , 0 W and 14"0 N from the center of the host 


galaxy NGC 4414 (Ciaba ttari et al.l 2013) at a distance 
of d ~ 16.6 ± 0.4 IVlpc (Freedman et al. 2001). The 
spectroscopic clas sification ot SiN 20 1 3dt be ing ot type lib 
was provided by |Ciabattari et al.| ([2013]) on 10.8 June 
2013. Based on the comp arison oFnght curve ev olution of 
SN 2013df and SN 1993J Van Dyk et al. (2014) estimate 
that the supernova occurred on JD 2,456,447.8 ± 0.5 or 
June 4.3. Throughout this article, we use this to be the 
date of explosion. 

We first observed SN 2013df with Jansky Very Large 
Array (VLA) on 2013 Jun 14.0 UT at the VLA C 
band. We diemnot detect the supernova on this occa¬ 
sion and also on our subsequent observation on 2013 
June 26.1 UT, about 10 and 22 days after the super¬ 
nova, respectively. On our third attempt on 2013 July 
5.9 we detected a bright radio source coincident with 
the optical position at a(J2000) = 12 h 26 m 29.33 s and 
J(J2000) = +31°13'38.3" (±0.1" in each coordinate) 

6 The Jansky Very Large Array is operated by the National 
Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 


TABLE 1 

VLA RADIO FLUX DENSITY MEASUREMENTS OF SN 2013DF a 


Date 

(UT) 

MJD 

Frequncy 

(GHz) 

F±a b 

(mJy) 

Array 

Config. 

2013 Jun 14.00 

56457.00 

4.8 

< 

0.063 

c 



7.1 

< 

0.051 


2013 Jun 26.07 

56469.07 

4.8 

< 

0.081 

c 



7.1 

< 

0.069 


2013 Jul 05.86 

56478.86 

8.6 

0.651 

± 

0.046 

c 



11.0 

0.732 

± 

0.034 




13.5 

0.647 

± 

0.041 




16.0 

0.615 

± 

0.033 


2013 Jul 21.09 

56494.09 

13.5 

2.107 

± 

0.030 

c 



16.0 

1.928 

± 

0.029 




19.2 

1.934 

± 

0.031 




24.5 

1.406 

± 

0.043 




30.0 

0.939 

± 

0.044 




43.7 

0.351 

± 

0.092 


2013 Aug 11.11 

56515.11 

5.0 

1.552 

± 

0.043 

c 



7.1 

2.639 

± 

0.041 




13.5 

2.356 

± 

0.035 




16.0 

2.145 

± 

0.035 


2014 Jan 27.36 

56684.36 

1.5 

1.531 

± 

0.100 

BnA 



5.0 

2.303 

± 

0.023 




7.1 

1.748 

± 

0.020 




13.5 

0.909 

± 

0.023 




16.0 

0.773 

± 

0.022 




33.5 

0.279 

± 

0.038 



a The quoted upper limits are 3cr 
b image rms 


with flux density of / y = 0.65 ± 0.05 mJy at 8.6 GHz. 
This, and subsequent observations from 1.5 to 44 GHz, 
are summarized in Table 1. 

All observations were taken in standard continuum ob¬ 
serving mode with a bandwidth of 16 x 64 MHz. During 
the reduction we split the data in two sub-bands (8 x 64 
MHz) of approximately 1 GHz each. We calibrated the 
flux density scale using observations of 3C 286, and for 
phase referencing we used the calibrator 4C 21.35 (JVAS 
J1224±2122). The data were reduced using standard 
packages within the Astronomical Image Processing Sys¬ 
tem (AIPS). 

2.2. Radio Spectral Energy Distribution of SN 2013df 

The overall spectral energy distribution (SED) of 
SN 2013df appears to be brightening with time at least 
until the observations of August 11, 2013 or within 
first 70 days after the SN. The peak spectral luminosity 
(d = 16.6 Mpc) L Vp = 3 x 10 26 erg/s/Hz on 11 August, 
2013 is about 4 times brighter than it was for the first 
detection only a month earlier on 2013 July 5.9. Interest¬ 
ingly, the spectral index f3 (for oc v?) is rather steep 
at early times (/3 = —1.6 ± 0.2 on 21st July 2013) and 
becomes shallower as it ages (/? = —0.9 ± 0.1 on 27th 
January 2014). 

This trend is unusual for the radio emission from SN. 
As the SN shock wave sweeps through the circumstellar 
medium (CSM) around a massive star, it shock acceler¬ 
ates the CSM electrons which generates the synchrotron 
emission. For the wind density profile (n(r) oc r~ 2 ) found 
around massive stars the total amount of matter swept 
up by the shock wave in a unit distance, or equivalently in 
unit time for a non-decelerating shock-wave, remains the 
same. As a result, the constant spectral peak brightness 
routinely observed in radio supernovae is the hallmark 
of the wind density profile. The evolving spectral index 
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is also suggestive of a radiation process in addition to 
the synchrotron emission and is evidently predominant 
at early times. These properties are considered in detail 
below. 


2.3. X-ray Observations using Swift and Chandra 
2.3.1. Swift-XRT 

The X-Ray Telescope ( XRT, IBurrows et al.||2 005D on¬ 
board the Swift satellite (Gehrels et al.' 2004|) started 
observing SN 2013df on 2013 June 13 (dt = 8 days since 
explosion, PI Routing). We analyzed the XRT data using 
HEASOFT (v6.15) and corresponding calibration files. 
Standard filtering and screening criteria were applied. A 
fading X-ray source is clearly detected at the position of 
SN 2013df that we associate with the SN shock interac¬ 
tion with the medium (Fig. [I]). 

Figure H] shows the temporal evolution of the X-ray 
emission from SN 2013df. We find a decay consistent 
with a power-law oc t~ a with a = 0.71 ± 0.16 (1 a c.l.). 
The contribution of diffuse X-ray emission from the host 
galaxy to the Swift-XRT data has been estimated and 
subtracted b y usin g the Chandra X-ray observations be¬ 
low (section 2.3.2) where the point-like emission from 
SN 2013df is well resolved. 

A combined spectrum made from the entire Swift- 
XRT data set is well modeled by an absorbed power- 
law with hard photon index T = 1.58 ± 0.06 (ler) and 
no evidence for neutral hydrogen absorption intrinsic to 
the host galaxy, in agreement with the results from the 
spectral analysis of Cha ndra o bservations which we dis¬ 
cuss in the next section 2.3. 2| The Galactic neutral hy¬ 
drogen column density in the d irectio n of SN 2013df is 
Nhmw ~ 1-60 x 10 20 cm~ 2 ( Kalberla et al.|[2005 ). 


2.3.2. Chandra 


We initiated deep X-ray follow up of SN 2013df with 
the Chandra X-ray Observatory on 2013 July 10, cor¬ 
responding to 5t = 36 days since explosion (PI Soder- 
berg). Data have been reduced with the Cl AO soft¬ 
ware package (version 4.6) and corresponding calibration 
files. Standard ACIS data filtering has been applied. In 
this observation which lasted for 9.9 ks we clearly de¬ 
tected SN 2013df, with significance > 50 er. As for the 
Swift-XRT observations, the spectrum is well modeled 
by an absorbed power-law with a hard photon index 
r = 1.26±0.25 (lcr). We again find no evidence for intrin¬ 
sic neutral hydrogen absorption, with a 3cr upper limit 
of Ni-f tnl < 2 x 10 21 cm -2 . The complete X-ray light- 
curve of SN 2013df comprising Swift-XRT and Chandra 
observations is presented in Figure [T] 

The hard X-ray spectrum can be fit by a thermal 
bremsstrahlung model with kT > 5 keV (Figure [2]). Due 
to relatively small number of counts and lack of coverage 
in hard X-ray band, we could establish only a lower limit 
on the plasma temperature. 

3. INVERSE COMPTON MODEL 


The main radiation process powering radio SNe is syn¬ 
chrotron emission. A shock wave is an efficient way to 
convert bulk kinetic energy into thermal motion of the 
shock accelerated particles. The SN shock wave accel¬ 
erates electrons to relativistic speeds which then gyrate 
in the shock amplified magnetic field and radiate syn¬ 
chrotron radiation. The electrons are usually assumed 


to have a power-law energy distribution, n e ( y e ) oc 7 “ p , 
with index p and Lorentz factor y e of individual electrons 
as a proxy of electron energy. This form of the electron 
distribution is evident from the synchrotron radio spec¬ 
trum. 

Two radiative cooling processes can affect this particle 
distribution: synchrotron cooling and inverse Compton 
scattering (IC). The synchrotron, as well as IC, power 
of an individual electron is strongly dependent on its 
Lorentz factor, P v oc y 2 . As a result, high energy elec¬ 
trons tend to cool down faster. The cooling timescale for 
electrons emi tting synch rotron rad iatio n at frequency v 
is given by Bjornsson & Fransson (2004) as 


t sc = 1.7 x 10 2 B 3 ^ 2 ^ 10 1 ^ 2 days (1) 

when viq = J//10 10 Hz FI Thus, the synchrotron cooling 
will be dominant only when t syn is the smallest timescale 
in the system i.e. t sc <C (f, U c ). 

The particle distribution could also be affected when 
high energy electrons cool down by scattering off optical 
photons due to photospheric emission from the SN. In 
the process, electrons cool down by imparting a fraction 
of their energy to the photons. The cooling time scale for 
this process depends on the reservoir of optical photons. 
For a SN having bolometric luminosity L and size r the 
average photon energy density would be given by U p h = 
L/ 47rr 2 c. Using this the cooling time scale is given by 


ni/2 r 2 

t ic = 1.5-jyf days (2) 

L 42 o 

There is another way of demonstrating the importance 
of IC as the dominant cooling process over synchrotron 
emission. It can be shown that the cooled down elec¬ 
trons have a steeper Lorentz factor distribution n e (y e ) oc 

7 e which is similar in the case of IC and synchrotron 
cooling, since both the processes have similar dependence 
on electron Lorentz factor: P u oc 7 1{Ub + U p h)- If the 
local magnetic field energy density ( Ub ) dominates over 
the local photon energy density ( U p h ), synchrotron cool¬ 
ing will be stronger than IC cooling and vice-versa. The 
effective electron distribution, and therefore the resultant 
radiation spectrum, will display only one break but its 
evolution with time will be determined by the dominant 
cooling process responsible for the break. 

Let "f sc and 7 j c be the characteristic Lorentz factors 
where the distribution of emitting particles steepens due 
to the synchrotron and IC cooling, respectively: 

7sc = 8.9x10 3 B- 2 f - 1 (3) 

Tic = 134 r 2 5 L 42 t d (4) 

The characteristic synchrotron frequency of emission due 
to an electron of Lorentz factor q e is v e = 1.2 x 10 6 f? 7 2 
Hz. Using this, the ratio of synchrotron and IC cooling 
break frequencies would be ( v sc /v ic ) oc L 2 2 B~ 4 rf 5 4 . We 
assume that the magnetic field energy density is a con¬ 
stant franction (e^) of the thermal energy density behind 
the shock and eb ~ 0.1. The thermal energy density 
behind the shock wave is given by U = ( 9/8)nm p /3 2 c 2 . 


1 Throughout this paper we will follow the notation X a = 
X/10 a , with quantities expressed in cgs units 
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Fig. 1.— Left panel : X-ray emission from SN 2013df is clearly detected by the Swift-XRT (0.3-10 keV). The image comprises observations 
collected starting June 13, 2013 until August 6, 2013, total of 49.5 ks (i.e. ~ 8 — 60 days since explosion). Right panel : The X-ray emission 
from SN 2013df as captured by the Swift-X RT and Chandra X-ray Observatory (CXO) shows a temporal decay consistent with a power-law 
t~ a with best-fitting a = 0.71 zb 0.16 (1 a c.l.). 



Frequency (10 15 Hz) 

Fig. 2.— Radio to X-ray SED of SN 2013df at the time of 
the Chandra observations. The X-ray flux is clearly in excess to 
the synchrotron model that best fits the radio observations (grey 
dashed line) and can be explained by thermal bremsstrahlung emis¬ 
sion with temperature T > 5 keV. Observations at E > 10 keV 
such as might be done by NuSTAR is necessary to constrain the 
temperature of the emitting plasma. 


Expressing the CSM density n in terms of mass-loss rate 
(M in units of 10 -6 M 0 /yr and wind velocity v w in units 
of lOOOkm/s) and r = /3ct for a freely expanding shock 
wave, gives the ratio of the break frequencies to be; 


— = 2.6 x ICE 


T 2 
^42 




1 C R ,-1 


M- 6 

Vw, 3 


-2 


(5) 


It is important to note that this ratio is independent 
of the source size. We can therefore state a general re¬ 
sult that for low to intermediate mass-loss rates (M < 
10 ““Mg/yr^ the dominant cooling process in the radio 
SNe would be inverse Compton rather than synchrotron 
cooling. 

The only time dependence in equation [5] is through 
the evolution of photospheric emission expressed as L 42 . 
The optical peak brightness is typically reached around 
20 days after the SN. Since v sc oc td, which follows from 
Equation [4] and subsequent discussion, it can be shown 
that the synchrotron cooling break is already out of the 


radio frequency bands by the time SN is at its peak bolo- 
rnetric luminosity. Therefore, the synchrotron cooling 
break would not be visible in typical radio SNe. 

Having demonstrated that the early steep radio SED 
is likely due to the IC scattering of optical photons, we 
now proceed to discuss model fitting and estimation of 
physical parameters for SN 2013df. 


4. MODEL FITS, SPECTRAL AND PHYSICAL 
PARAMETERS 

As discussed in the previous section IC cooling modi¬ 
fies the electron distribution and consequently the syn¬ 
chrotron emission. Evolution of this electron distribu¬ 
tion with time depends on the bolometric luminosity as 
well as the size of the SN. In order to parametrize this 
evolution we have characterized the SED using the syn¬ 
chrotron self-absorption break frequency (v a ), the char¬ 
acteristic synchrotron frequency due to electrons with the 
minimum Lorentz factor ( v m ) and the inverse Compton 
cooling break frequency (za c ), as well as the spectral peak 
f Va of the synchrotron spectrum. In the appendix, we 
provide spectral and temporal scalings for the evolution 
of these spectral breaks. 

To fit all the radio observati ons co ll ective ly we fo llowed 


a procedure similar to that of Kamble et al. (20141 and 
scalings from the Appendix here. Ey using least square 
minimization we estimated values for the spectral peak 
f Va , break frequencies v a ,Vi C , frequency of free-free ab¬ 
sorption vs where optical depth to free-free absorption 
Tff = 1 (for definition see Equation [h] below), electron 
distribution index p and deceleration parameter m. The 
bolometric luminosity of SN 2013df was used as an input 
to determine the evolution of v ic . 

The statistical uncertainties in flux densities measured 
and listed in Table [T] are the respective image RMS. Vari¬ 
ations in the local weather at the telescope can introduce 
additional errors to the measured image RMSs. Fur¬ 
thermore, calibration of the flux and phase calibrators 
can add further systematic errors of a few percentage. 
Therefore, during the model fitting, about 15% of the 
total error was added in quadrature to account for these 
sources of errors. 

The radio SED indicate that the frequency v m < v a 
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Fig. 3.— Radio SED of SN 2013df at different epochs. The lines indicate the best fit models: SSA as the dotted lines, FFA (g> SSA as 
the dashed lines and IC (g> FFA ® SSA as the solid lines. 


during all the observations and in any case is below the 
observing band. As a result, we can constrain the spec¬ 
tral break i/ m <Cl GHz. 

The best fit value of spectral peak flux density is de¬ 
termined to be f Va = 48 ± 5 mJy at the self-absorption 
frequency v a = 3.3 ± 0.5 GHz on 21st July 2014 or at 
t = 46.8 days after the burst and the best fit reduced 
y 2 = 5.2 per degrees of freedom (DOF) for 14 DOF. 
We note that a significant contribution to the y 2 comes 
from the earliest epoch of radio detections on 05 July 
2013. Dropping those observations from the fit results 
in a significant improvement with reduced y 2 = 1.7 for 
10 DOF. The resultant model SEDs are compared with 
observations in Figurepjand the light-curves in Figure |4j 

The Figures [3] and [4] highlight the importance of man 
vidual emission and suppression processes such as SSA 
(synchrotron self-absorption), FFA (free-free absorption) 
and IC (inverse Compton). At early epochs, the model 
synchrotron emission is under-estimated and can not an¬ 
ticipate the subsequent rise very well. The rising SEDs 
suggest an a ccelerating shock wave with a r = 1.12 (see 
Equatior jA2[ ). The best fit parameters obtained using 
the SSA model alone (reduced y 2 = 3.8 for 16 DOF), 


f Va = 1.8 mJy and v a = 15.6 GHz result in an unre¬ 
alistically slow shock wave v s h = 0.02 c. This is signifi¬ 
cantly slower than the velocities estimated by the emis¬ 
sion lines in optical spectra and typical ejecta velocities 
of v e j ss 10,000 km/s in SNe. This, therefore, clearly 
indicates that the free-free absorption must be playing 
an important role in SN 2013df. 

Inclusion of free-free absorption improved the esti¬ 
mates of the shock velocity consistent with the typical 
values, but not necessarily the quality of the fits (reduced 
y 2 = 6.8 for 15 DOF). The SSA and SSAoFFA models 
however result in shallower SEDs which do not agree well 
with observations at early epochs, most notably on 21st 
July 2013. The steep spectral shape should be related 
to the underlying steep electron energy distribution and 
therefore could be the result of a cooling process. As 
discussed in section § [3] inverse Compton cooling is likely 
to play the significant role. We therefore believe that 
the steep SED is the result of inverse Compton cooling 
process. 

At late epochs, the radio emission at GHz frequen¬ 
cies is expected to be less affected by the suppression ef¬ 
fects: a) the synchrotron self absorption (v a ) falls below 
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Time (MJD - 56447.30) 


Fig. 4. — Radio light curves of SN 2013df at various radio-frequency bands. The lines indicate the best fit models: SSA as the dotted 
lines, FFA 0 SSA as the dashed lines and IC 0 FFA 0 SSA as the solid lines. 


GHz bands, b) effects of free-free absorption diminishes 
rapidly with increasing SN size as the optical depth to 
free-free absorption Tq oc r -3 and c) inverse Compton 
cooling becomes insignificant as the reservoir of optical 
photons is drained off with diminishing bolometric lumi¬ 
nosity. As a result, all three model estimates roughly 
agree with each other on the last epoch of observations 
on 27th January 2014 i.e. about 237 days after the SN. 

For the electron energy distribution index, the fit con¬ 
verged to p = 2.7 ± 0.3. The observations before 50 
days show significant effects of inverse Compton cooling 
through the steep spectral index jz - ' 3 where /3 = 1.6±0.2. 
As discussed in the appendix this spectral regime is re¬ 
lated to the electron distribution index through /3 = p/2 
and is consistent with the measured p = 2.7 ± 0.3. At 
later epochs, when the inverse Compton subsides, opti¬ 
cally thin synchrotron spectrum is related to the electron 
distribution as /3 = (p— 1)/2. The observed /3 = 0.9±0.1 
is therefore consistent with it. 

The shock wave expansion is best f it as r = ro(t/to) ar 
with a r = 0.84 ± 0.1 (see Equation A2|. This is con¬ 
sistent with the expectations for a compact, radiative 
envelope star where the outer density profile should tend 


to a power law i n radius, p m cx r 30,2 gi ving a r = 0.88 


Chevalier| 

1982 

Matzner & McKeeJ1999 

). Our estimate 

a r = (J.84 J 

= 0.1 is also consistent with a convective outer 


envelope ( p sn oc r 12 ) where one expects a r = 0.9. 

4.1. X-rays as Thermal Bremsstrahlung emission 

X-ray emission in SNe could be produced by one of the 
three radiation processes: synchrotron, inverse Compton 
and bremsstrahlung. Using the radio spectral index of 
/3 = —1.6 based on observations and modeling of the 
radio emission, we expect the synchrotron emission in 
X-ray bands to be more than two orders of magnitude 
fainter as shown in Figure [2] 

Inverse Compton scattering of the photospheric opti¬ 
cal emission in to the X-ray bands by relativistic non- 
thermal electrons behind the forward shock is another 
possibility. The X-ray spectrum produced by such scat¬ 
tering would be soft. The spectral photon index of the 
observed X-ray emission in SN 2013df (T = 1.3) is, how¬ 
ever, significantly hard and is incompatible with a spec¬ 
trum dominated by Comptonization. Instead, such a 
hard spectrum is consistent with the expectations from 
thermal bremsstrahlung emission. 
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SN 1993J was another type lib supernova with bright 
X-ray emission detected for up to sever al tens of days. 


Fransson, Lundqvist, & Chevalier (1996) modeled it us¬ 
ing the thermal bremsstrahiung emission from the for¬ 
ward and reverse shocks that propagate into the stellar 
wind and ejecta, respectively . Fo l lowing a d etaile d ana- 
lytic argument Fransson, Lundqvist, & Chevalier (1996) 
showed that the reverse shock was radiative and most of 
the X-ray emission during the first few months originated 
in the forward shock. They estimated the pre-SN mass- 
loss rate to be ~ 4 x 10 -5 M 0 /yr for a wind velocity of 
10 km/s. 

Following the p roce dure outlined in Fransson/1 
Lundqvist, & Chevalier ([1996), we come to the similar 
conclusions for SN 2013df as those for SN 1993J. In Fig¬ 
ure [5] we show a range of mass-loss rates and forward 
shock velocities to compare the resultant X-ray flux due 
to thermal bremsstrahiung emission with observations. 

As pointed out by |Fransson, Lundqvist, fe Chevalier| 
(1996) the X-ray emission may originate either in the 
forward shock that is sweeping the dense pre-SN wind 
or in the reverse shock that forms at the contact dis¬ 
continuity and propagates into the denser ejecta. For 
homologous expansion, as is approximately the case for 
SNe, it is convenient to express the ejecta density profile 
is as a function of velocity, v, as p sn oc v~ n . The reverse 
shock te mperatu r e may depen d on thi s density pr ofile. 
Following Fransson, Lundqvist, & Chevalier (1996), we 
find that the temperature of the reverse shock T rs > 10 7 
K and that of the forward shock Tf s > 10 9 K for reason¬ 
able shock velocities v s h > 10 4 km/s and density profile 
of the ejecta n < 15. 

If the reverse shock is radiative then a cool shell forms 
between between the reverse and forward shocked ma¬ 
terial which absorbs all the X-rays from the reverse 
shock. For T > 2 x 10' K, cooling is dominated by 
bremsstrahiung emission and below that it is dominated 
by line emission. At such temperatures reverse shock 
remains radiative for a significant amount of time. For 
M > 8M 0 /yr, n > 12 and v s h < 21,000 km/s we find 
that the reverse shock remains radiative for at least 50 
days after the SN which spans the entire duration of our 
X-ray observations. Thus, all X-rays from the reverse 
shock are absorbed and most of the observed X-ray flux 
is due to the forward shocked material. 

Chandra observation provides the best constraints on 
the model parameters. The red shaded area in Figure 
[H] corresponds to 3cr contour bounded by (6.1 < M < 
10.1) x lO~ 5 M 0 /yr (v w = 10 km/s) with the forward 
shock velocity being (0.05 < /? < 0.1). The correspond¬ 
ing mass-loss rates for the 5(7 contours shown in larger 
shaded area are (4 < M < 11.1) x 10 -5 M 0 /yr, v w = 10 
km/s for similar shock velocities. The solid red line 
corresponding to the best model estimate for reproduc¬ 
ing Chandra observations gives M = 8.3 x lO“ 5 M 0 /yr, 
v w = 10 km/s and /3 = 0.07. 


4.2. Model Parameters from Radio Emission 

Another effect of the stellar mass loss is in the form 
of absorption of radio emission by the ionized wind sur¬ 
rounding the SN. The free-free absorption optical depth 



Time [days] 


Fig. 5.— Thermal bremsstrahiung emission in X-rays. The 
Chandra and Swift observation are shown as blue data points. 
The modelled thermal bremsstrahiung emission is shown by the 
shaded region. The lower bound corresponds to the shock velocity 
P = 0.05 c, mass-loss rate parameter A * = 4 and the upper bound 
to ft = 0.1, A* = 11. 


of this unshocked gas at frequency v is given by 


Tff 


1.9 A 9 
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where /3 is the shock velocity in units of the speed of light 
and and A* is the mass loss rate in terms of a reference 
value of 5 x 10 13 gm/cm. The mass loss parameter A = 
Af/47rI4jgm cm -1 and our reference value of 5 x 10 13 
gm/cm is attained for M = lO” 5 M 0 /yr and v w = 10 
km/s. 

At t = 46.8 day after the SN, our best fit model gives 
Tff = 1 at v = 10 A 1 GHz. Equation [G] can then be used 
to get the mass loss rate parameter in terms of shock 
velocity (3 a s A* = 584 d 3 / 2 . 


-3 


Following Kamble et al. (2014), we express the syn¬ 
chrotron frequency and flux m terms of unknown model 
parameters: 
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where e e and es are the fractions of thermal energy be¬ 
hind the shock that goes into accelerating electrons and 
magnetic fields, respectively. 

Using the best fit spectral parameters determined from 
radio emission modeling described in the previous sec¬ 
tion one can simultaneously solve Equations [8] and [6] for 
physical parameters (3,e e and cb- Using this approach 
we estimate the shock velocity to be /3 = 0.07 ± 0.01 
which is faster that the velocity of the fastest ejecta esti- 
mated from the optical line-widths of about 10,000 km/s 


(Morales-Garoffolo et al.|2014| Ben-Ami et al. |2014[ Van 
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L)yk et al.|2U14 Maeda et ai.|2015). This estimate is also 
consistent with the thermal bremsstrahiung emission ob- 
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served in X-rays (see Figure[l]). For SN 1993J the velocity 
estimated from optical lines was about 19, 000 km/s, and 
a shoc k velocity was est imat ed to be 21,600 — 24,700 
km/s (Fransson & Bjornsson 1998b in agreement with 
the VLB! measurements ot Bartel et ai.j (1994 2002), 
which gives /3 = 0.073, same as that estimated above 
for SN 2013df. The SN mass-loss rate is then estimated 
to be M = 10.5 ± 3.0MQ/yr for the wind velocity of 
v w = 10 km/s, in agreement with the estimates from the 
X-rays. It is reassu ring that in an independent approach 
Maeda et al. (2015) reach roughly similar conclusions us¬ 
ing late time optical evolution of SN 2013df. 

Similarly, we determine e e = 0.2 ± 0.02 and cb = 
(1.0 ± 0.1) x 10 -3 . While this determined value of e e is 
similar to those generally found in SNe, we find a signif¬ 
icant deviation from energy equipartition between mag- 
netic fiel d s an d acce ler ated particles. For comparison, 
Fransson & Bjornsson ((1998) deduced es = 0.14 for SN 
1993J. The high mass-loss rate combined with low cb 
may be the reasons for unusual radio properties of SN 
2013df. The only other SN for which inverse Compton 
effect was shown to be of importance in shaping radio 
and X-ray emission was SN 2002ap, a type Ic supernova. 
It is interesting to note that a low value of e# = 2 x 10~ 3 , 
similar to the pr esent case of SN 2013df, was deduced for 
SN2002ap also (Bjornsson & F ranssonf200 4). 


4.3. Inverse Compton Cooling Break in Radio 

The break in the energy distribution of relativistic elec¬ 
trons due to the IC cooling corresponds to a spectral 
break vie discussed in section [3j Following that discus¬ 
sion, it can be shown in terms of the physical parameters 


v IC = 1.4 x 10 10 ^ 

-^42 


M. 


-5 


Vw, 1 


t d y/e B -1 Hz (9) 


Due to its strong dependence on shock velocity, vie 
favors a slow shock wave and lower mass-loss rates. We 
obtained vie = 0.6 GHz at to = 46.8 days for the best 
fit to the observations. The bolometric luminosity of 
SN 2 013d f at that epoch is sa (0.8 ± 0.3) x 10 42 erg/s 
(Morales-Garoffolo et al. 2014). For these observables 


one can estimate a combination /3e 2 J 19 = 0.01. For the 
fastest line velocities /3 = 0.04 which gives es = 10” 5 . 
However, due to the strong dependence of vie on other 
quantities, especially Lt> 0 i and /3, this comes with a large 
uncertainty es = (1 ± 2) x 10 -5 . 


4.4. Inverse Compton emission in X-rays 

The relativistic electrons which generate the syn¬ 
chrotron radiation and powers the radio emission in SN 
shock wave, will naturally up-scatter the optical photons 
that emerge from the photosphere of the SN. If the seed 
photons or the relativistic electrons have a wide energy 
distribution, the spectrum of the up scattered inverse 
Compton emission will also be wide. One can, how¬ 
ever, estimate the total contribution of IC component in 
a quasi-monochromatic X-ray band. The resultant X-ray 
emission will be related to the energy distribution of elec¬ 
trons. The non-thermal electron distribution will pro¬ 
duce non-thermal X-rays. In the process of IC scatter¬ 
ing, the electrons lose their energy. Higher the electron 
Lorentz factor, faster it loses energy to inverse Compton. 


This, in turn, changes the shape of the electron distribu¬ 
tion at higher Lorentz factor, making it steeper than the 
electron distribution injected at the shock front. 

Radio observations of SNe are an excellent probe of the 
underlying electron distribution and magnetic field be¬ 
hind the shock. One can construct bolometric lightcurve 
of SNe using optical observations. Thus, between ra¬ 
dio and optical, one has all the necessary information 
required to understand and estimate the inverse Comp¬ 
ton emission. Below we use these observations to self- 
consistently model the IC emission. 

From the previous section, radio observations on 21st 
July 2013, or at supernova age of to = 46.8 days, in¬ 
dicate B « 0.8 G. Therefore, the magnetic field en¬ 
ergy density behind the shock is Ub = 0.02erg/cm 3 . 
Radio observations also indicate the size of the shock 
wave to be r = 10 16 cm. From the optical observa¬ 
tions we estimate the bolometric luminosity of the SN 
to be Lb 0 i ss 0.8 x 10 42 erg/s on the the same epoch. 
Thus, energy density in seed photons turns out to be 
U p h ~ 0.03 er g /c m 3 . Furthermore, £(p = 2.7) = 3.06 
(see Equation B1 in the appendix), and f Va = 48 mJy at 
the self-absorption freq uency v a = 3.3 GHz. Using these 
values in Equation |B1| and p = 2.7 as determined from 
radio observations, we determine F v {\ keV) = 2.3 x 10 -6 
mJy. This is more than an order of magnitude fainter 
than the observed X-ray emi ssion . 

Similarly, using Equation |B1| and p = 2.7, one can 
estimate the time evolution of IC X-ray flux to be 
F„(lkeV) oc L bol (t) 0 963 x t -0 925 . The expected light 
curve of the IC X-ray emission in 1-10 keV band is sig¬ 
nificantly fainter. Overall, the contribution of the inverse 
Compton emission to observed X-ray emission is negligi¬ 
ble. 

In other words, the observed X-ray emission is too 
bright at all epochs than what is expected from the in¬ 
verse Compton emission. One can therefore conclude 
that the X-rays from SN 2013df must be due to another 
strong emission component. X-ray brightness estimates 
due to thermal bremsstrahlung, as shown in the previous 
sections, appear to be reasonable. 


5. DISCUSSION 

5.1. Radio SNe - a global view 

It is instructive to view the SN 2013df in comparison 
with radio SNe in general. The peak Luminosity v/s 
peak time plot offers the most informative way to project 
the radio SNe. The positions of radio SNe in such a plot 
are shown in Figure |6( which i s an update of the previous 
plots in Chevalier (1998 );|Chevalier, Fransson, & Nymark 
(2006). Different SN types have been color coded tor the 
ease of distinction - Ib/c: blue pentagons; lib: red circles; 
IIP: magenta squares; Iln: green diamonds; IIL: cyan 
triangles and II (SN 1978K, 1981K & 1982aa): yellow 
octagons. 

While both the co-ordinates of SNe on this plot are 
purely observables, the underlying lines depend, though 
weakly, on certain model assumptions. The blue dotted 
lines of constant shock velocity assume energy equipar¬ 
tition between synchrotron emitting electrons and mag¬ 
netic field, and that there is no strong cooling of the elec¬ 
trons. The lines also assume a uniform electron energy 
distribution index of p = 3. Despite these assumptions 
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Fig. 6. — Peak spectral radio luminosity v/s product of the peak time and frequency for radio SNe of various types - Ib/c: blue pentagons; 
lib: red circles; IIP: magenta squares; Iln: green diamonds; IIL: cyan triangles and II (SN 1978K, 1981K &; 1982aa): yellow octagons. The 
similarity between SN 1993J and SN 2013df is clearly evident in this plot. 


the plot provides an easy way of comparison between dif¬ 
ferent SNe as well as insights into different classes of SNe 
and their progenitors. 

The slowest to fastest moving objects appear from right 
to the left of the plot. The blue dotted lines indicate 
mean velocities of the radio shocks if SSA is the respon¬ 
sible for the spectral peak. The relativistic synchrotron 
sources such as gamma ray bursts appear to the far left 
of the plot in the dark green shaded area. The relativistic 
SNe, such as SN 2009bb, appear in the light green shaded 
area labelled as the engine driven SNe. As one moves to 
the right of the plot synchrotron emission is increasingly 
suppressed by the FFA in the progenitor wind. 

The effect of FFA is often noticeable in radio observa¬ 
tions where the SED falls faster than in the SSA towards 
low radio frequencies. Consequently, the rise of radio 
light curve is also faster at early times in the presence of 
FFA. Our early non-detections of SN 2013df at C band 
clearly bring out the presence of FFA. Ignoring FFA re¬ 
sults in the estimates of reduced shock-velocity, which 
often fall short of photospheric velocities measured from 
optical line emission. Therefore, we use nominal a veloc¬ 


ity of < 10,000 km/s as a distinction between FFA and 
SSA region in the luminosity plot (Figure[6]). 

5.2. Progenitors of SNe lib - extended versus compact 

Figure [6] captures the essence of radio emission in SNe 
and portrays the diversity among various types of SNe. 
Radio SNe show a wide luminosity distribution spread 
over more than four orders of magnitude. Luminosity 
distribution of lib SNe is significantly narrow, confined 
within two orders of magnitude. But their distribution of 
peak times is relatively as wide as those of SNe Ibc. The 
distribution merges with those of SN Ib/c. As one moves 
to the right, the effect of FFA becomes increasingly no- 
ticeable in SNe li b as has be en clearly seen in SN 1993J 
( Fransson fc Bjorns son|1998), 20 13df (this work), 2011hs 
(Bufano et al.|2U14|) and 2U10P (iRomero-Canizales et al. 

It is, therefore, plausible that the broad distribution in 
rise times of SN lib may be an outcome of different in¬ 
trinsic properties of the progenitors or the circumstances 
surroundin g t h e SN. Indee d, such a propos al has been put 
forward by Chevalier & Soderberg (2010) based on dis- 
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tinct emission properties observed in radio, optical and 
X-rays. 

One clue to this might come from early rapidly evolv¬ 
ing optical emission often termed as the ‘cooling envelope 
emission’. Current understanding of this emission is that 
following the shock break-out from stellar surface, outer 
envelope of the star expands and cools. As the expanding 
photosphere reaches this stellar envelope the radiation 
stored in the envelope manages to escape. This ‘hot’ radi¬ 
ation results in the rapidly evolving UV/optical emission 
at early times lasting from a few hours to several days. 
The cases where such emission has been detected are 
only a few consisting of SNe of type lib: 1993J, 201 ldh, 
201 lhs, 201 lfu and most recently 2013df, and type Ic: 
SN 2006aj and a few cases of SN of type IIP. 

This however is not a complete picture. Observations 
are maturing rapidly in discovering SNe early, in carrying 
out rapid follow up and in UV capabilities. Rapid follow 
up reveals that SNe does not always develop such double 
peaked early optical emission. Despite observations at 
very early times, several SNe have been detected which 
do not show early optical emission to very deep limits 
which might be problematic to the theories of cooling 
envelope emission, especially if the SNe are expected to 
have similar progenitors. It may, therefore, be insightful 
to tie up such properties with multiband long term follow 
up observations. 

A simple but surprising trend seem to emerge between 
this early optical emission and slowly evolving SNe. And 
at the moment we concern ourselves with of those of type 
lib only. All SNe in the FFA region of Figure [G] exhibit 
this early optical emission whereas others do not ! As 
mentioned before, SNe 1993J, 2011hs, 2013df are the 
prime examples of this optical emission, they all show 
strong FFA absorption in the radio emission and reach 
radio peak brightness comparatively late (> 100 days). 
SNe 2010P which also falls in the same FFA category 
was never observed in UV/optical so we can not com¬ 
ment on that. Opposite of this is SN 2011fu, which ex¬ 
hibited the early optical peak but was never followed up 
in radio so it is not possible to comment if it also had a 
strong FFA component. On the other hand, SNe without 
strong FFA do not show this optical emission to deep lim¬ 
its despite early observations (SNe 1996cb, 2008ax and 
201 lei). Early UV/optical observations have not been 
reported for SN 2001gd and 2001ig. SN 201 ldh is an in¬ 
teresting intermediate case, where the early optical emis¬ 
sion was detected in one band and was significantly short 
lived. Also, its contribution to the bolometric luminosity 
is negligible. The progenitor of this SN has however been 
identified and we discuss that issue later. 

These findings are put in persecutive and more quan¬ 
titatively in Figure [7] where we plot the product of the 
early peak bolometric luminosity and its epoch against 
the epoch of radio peak brightness. The right panel of 
Figure [7] shows X-ray emission in SNe lib. Interestingly, 
the bright X-ray emission correlates with radio proper¬ 
ties such as FFA, delayed radio peak, and double peaked 
optical emission. Rapidly evolving radio SNe are also the 
ones which display faint or none X-ray emission. Unlike 
SN 1993J and 2013df, X-ray emission of SN 201 ldh was 
non-thermal in nature and more than an order of mag¬ 
nitude fainter. 

Early UV/Optical observations of SNe provide inter¬ 


esting diagnostics about the SN progenitor, especially 
in cases where pre-SN images of the progenitor do not 
exist. But before that the physical theories need to be 
cross checked and calibrated against the cases where the 
progenitor has been identified. Undoubtedly, such de¬ 
tections will improve in future and we will have more 
candidates to test theories against. At the moment we 
have four SNe with identified progenitors, properties of 
which we discuss in the next section. 


5.3. SNe lib with identified progenitors 

Observations are not only getting quick in respond¬ 
ing to the SNe fast enough to catch the early emission 
but also probing deeper to match and identify the pro¬ 
genitors from archival images of nearby SNe. As a re¬ 
sult progenitors of four SNe of type lib have been iden¬ 
tified in pre-SN images. A YSG of radius ~ 600 R G 
was identified as the progenitor of SN 1993J. A simi¬ 
larly large YSG star has been identified as a progenitor 
of SN 2013df. A less extended, a YSG nevertheless, with 
R ss 2OO1?0 was found to have disappeared from the po¬ 
sition of SN 201 ldh and has therefore been concluded as 
the progenitor of SN 2011dh. With the surge in direct 
identification of SN progenitors it is now becoming in¬ 
creasingly possible to test some of the predictions about 
early emission properties of similar SNe. 

The question then is why should similar kind of pro¬ 
genitors have different emission properties such as in the 
cooling envelope which emerge from close to the stellar 
surface ? As pointed out before not every type lib SN 
shows cooling envelope emission, and there appears to be 
correlation with the progenitor size. Indeed, SN 2008ax 
did not display the cooling envelope emission. It was dif¬ 
ficult to fit either a single or binary supergiant or WR + 
supergiant system self-consistently as progenitors to the 
observed colors from the pre-SN images (Li et al. 2008; 
Crockett et al. 2008). The possible initial mass range for 
the progenitor turns out to be large, ss 10 — 28 M 0 . 

More progress is required on theoretical front in figur¬ 
ing out all the elements for the CEE theory to b e suc cess- 
ful. It has been pointed out by Nakar & Piro (20141 that 
standard stellar density profile does not lead to CE emis¬ 
sion. Instead, extended progenitors with non-standard 
density profiles are necessary. Such modifications of den¬ 
sity profiles in the extended envelopes may occur only in 
a subset of progenitors due to binary interaction, non¬ 
standard mass-loss, episodic ejections of outer envelopes 
and subsequent evolution etc. There are, however, no 
direct means of observationally confirming if such mod¬ 
ified density profiles actually exist elsewhere and more 
so under a given situation. Secondly, current theories 
do not take into account de tailed d yna mics or rad iatio n 
transport as pointed o ut by Rabinak fe Wax man (2011); 
Nakar & Sari (20101; Makar & Piro 12014 1. Different 
kinds of progenitor s, such as WR , BSC a na RSG, have 
been c onsidered by Nakar & Sari (2010); Nakar & Piro 
(20141 to estimate CE emission properties! One ditter- 
ence that is expected, and to which calculations from 
different groups agree, is that larger progenitor results in 
longer lasting CE emission. 

A credible theory must explain all the apparent con¬ 
tradictions and complexities as pointed out above. This 
raises an important question: is the double peaked op¬ 
tical light curve really due to the CE emission ? Alter- 


























Radio and X-ray Supernova 2013df 


11 


> 

D 
5 n 

c3 

O 

a 

o 

i_i 

B> 

(U 

a 

H 

X 

a 


io 47 





• 

96cb 

▼ 

• 

llhs 

• 

93J 

13df 

llei 




T 

- 08ax 



- 

T T 

lldh 





SN2013df 






A ^ 




<> 


20 


40 


60 


80 


100 120 


140 


160 


Epoch of Peak Radio Luminosity (days) 


100 

Days since Explosion 


Fig. 7.— Early optical emission and radio properties: A product of the early peak bolometric luminosity and its epoch are plotted 
against the epoch of radio peak brightness (Left panel). The quantity on the y-axis is a measure of the amount of energy a SN deposits in 
the envelope. The red points are the SNe with detected early optical emission. The blue upper limits are the SNe where early UV/optical 
observations have been reported without similar long lasting optical emission. The detections and non-detections of early optical peak are 
clearly separated signaling apparent dichotomy. X-ray emission in type lib: Compared to a sample of type lib SNe, SN 2013df shows 
luminous X-ray emission with strong similarities to SN 1993J, both in terms of luminosity and temporal decay rate (Right panel). Th e SNe 
with bri ght X-ray em ission are also the ones w hich show F FA and delayed peak brightness in radio. Ref erences: SN 1993J (Chandra et| 


|al.|2009|), S N 2001gd llPerez- Tor res et al.l20 05t, SN 2003be Soderb erg. Chevalier. Kulkarni, fe Fralll2006ll. SN 2004C (Dittman m prep.), 
SN 2UUSax jRoming et al.|2UU9|l, SN 2Ulldh jboderberg et a .|2U1Z||. 


native ideas should be explored in order to explain ob¬ 
servations satisfactorily. Interaction of shock-wave with 
binary companion is expected to leave some signature 
on the optical light curves. Such an idea has been ex¬ 
plored for SN la where a white-dwarf accretes matter 
from the companion ultimately leading to the SN. As 
the SN shock-wave interacts with the envelope of the 
companion dissipating kinetic energy and heating it up, 
it results in the emission that is expected to peak in UV 
bands on timescales of a few days which is roughly similar 
to that observed in some type lib. It is of considerable 
interest that a binary companion has also been either 
confirmed (SN 1993J, SN 201 ldh) or speculated to the 
progenitors of these SNe. 

SNe of type lib thus provide the best opportunities to 
study variety of phenomena: double peaked bolometric 
light curves and bright radio and X-ray emission. Due 
to the relatively higher rate of type lib compared to Ibc 
they are also best suited for direct detection of progeni¬ 
tors. With the possibility of such rich datasets they are 
ideal for tying together various observed properties and 
in building a detailed understanding. Correlating pro¬ 
genitors with radio properties may be more subtle than 
the simple dichotomy of compact and extended progen¬ 
itors and may be made more complex by the processes 
that affect stellar evolution such as stellar rotation, mag¬ 
netic fi elds, b inar y interactions and evolution in the late 
stages (Smith & Tombleson 2015). A more complete pic¬ 
ture can be built by mufti band observations and direct 
detection of progenitors in the pre-SN images. 


6. SUMMARY AND CONCLUSIONS 

We presented extensive radio and X-ray observations 
of SN 2013df in the context of type lib SN explosions. We 
also presented a model of inverse Compton cooling and 
estimated its effects on the evolution of radio emission 
from SN 2013df. Our main results can be summarized as 
follows: 


• With bolometric luminosities within a factor of 
a few of each other, and similarly extended lu¬ 
minous cooling envelope emission and compara¬ 
ble mass-loss rate estimates from the radio and X- 
ray emissions, SN 2013df appears similar to that of 
SN 1993J, in several ways. This close association is 
portrayed most easily in Figure [6 Indeed, the pro¬ 
genitor of SN 2013df appears tohave be en similar 
to that of SN 1993J in mass and radius (Van Dyk 
et al.|[2014[ ). 


• SN 2013df showed a peculiar soft-to-hard tempo¬ 
ral evolution of the optically thin part of the radio 
spectrum. We interpret this behavior as free-free 
absorbed synchrotron radiation where the under¬ 
lying electron distribution responsible for the ra¬ 
dio emission is modified by Inverse Compton cool¬ 
ing at early times (ref to Fig. 3). This model 
provides a reasonably good match to the obser¬ 
vations for t > 40 days and leads to the picture 
of a shock wave propagating with modest velocity 
v s h = 0.07c into a dense environment, previously 
shaped by substantial mass loss from the progen¬ 
itor star M = 8 X lO -5 M 0 /yr (for wind veloc¬ 
ity Vw = lOkms -1 ). Comparable mas s-loss ra te 
was inf erred for the type lib SN 1993J (Fransson, 
Lundqvist, & Chev aher||199 6). 


• Self-consistent modeling of the radio emission of 
SN 2013df that includes synchrotron emission and 
self absorption, free-free absorption and inverse 
Compton cooling effects shows significant devia¬ 
tion from the generally assumed equipartition of 
shock energy between magnetic fields and acceler¬ 
ated electrons. For SN 2013df we infer e e ~ 200es. 


• In close similarity to SN 1993J, the X-ray emission 
from SN 2013df (Fig. 7) is in clear excess with re¬ 
spect to the synchrotron model that best fits the ra- 
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clio emission (ref Fig. 2). The hard X-ray spectrum 
is dominated by thermal bremsstrahlung emission 
with T > 5 keV and indicates a mass-loss rate 
(4 < M < 11) x 10~ 5 M(7)/yr, consistent with our 
results above. 

• Our modeling of the radio and X-ray radiation from 
SN 2013df further demonstrates the importance of 
the IC cooling process in shaping the electron dis¬ 
tribution responsible for the detected radio emis¬ 
sion in SNe at early times. 

• Finally, we observationally link the presence of 
bright cooling envelope emission in the optical/UV 
bands at very early times in type lib SNe to the 
later appearance of heavily free-free absorbed radio 
emission (Fig. 7), coupled with hard X-ray radia¬ 
tion from thermal bremsstrahlung (and hence large 
mass-loss rate of the progenitor star before explod¬ 
ing). The dichotomy of such emission properties 
displayed by type lib SNe suggest diversity among 
progenitors. 

• The current small number of type lib SNe well ob¬ 
served in the optical/UV, radio and X-ray bands 
both at early and at late times (Fig. 7) does 
not allow us to quantitatively test our conjectures. 
Clearly, a significantly larger sample of type lib SN 
explosions with extended multi-wavelength cover¬ 
age is needed. The case of SN 2013df demonstrate 
that multi-band observations provide crucial diag¬ 
nostics about the nature of SN progenitors. SNe 
of type lib are advantageous to carry out multi¬ 
faceted investigations including detections of pro¬ 
genitors in pre-SN observations, double peaked op¬ 
tical emission and non-thermal emission. 
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APPENDIX 


SYNCHROTRON EMISSION FROM INVERSE COMPTON COOLED ELECTRONS IN YOUNG SUPERNOVAE 


Synchrotron spectrum produced by the relativistic electrons behind a non-relativistic shockwave expanding into the 
surround medium of density profile n(r) oc r~ 2 could b e described as a bro ken power law with two char acteristic 
frequencies (v m , v a ) and the spectral peak F m . Following Soderberg, Chevalier, Kulkarni, & Frail ([2006) and Kamble 
et al. (2014) we express the size and speed of the shockwave evolving in time as 


r = r 0 x ( t/t 0 ) ar 
f3 = l3 0 x (t/to) ar_1 


(Al) 

(A2) 


where the shock speed is /3c = dr/dt. 


Spectral shape, Break Frequencies and Radio Light curves 

We identify the spectral break due to inverse Compton cooling of electrons as a characteristic frequency (v ic ). The 
evolution of this break frequency depends not only on the SN size but also on its bolometric luminosity, oc t aL . 
Below we describe two spectral regimes that are most relevant for the evolution of radio SNe. 


Case 1: v m <^v a <^ v ic 

In Case 1, the characteristic synchrotron frequency is below the self-absorption frequency, v m <C v a and therefore 
the synchrotron emission is suppressed due to the self-absorption. As a result, the spectral peak occurs at v p ss v a and 
the scalings simplify to 


_ (5p —2)—2a; r (2p —3) 

VpK, V a <X t P+ 4 , 


fv v « Av OC t 


l ■| \ 12p — 7 


Vir OC t 


Acx r —otL—3 


The temporal and frequency dependence of flux /„ at observed frequency v is then generalized by 

( ^ 5 / 2 t 2ar + 1 / 2 U m < V < v a 

f v OC i l /-(p- 1 )/ 2 t -[(5p-3)/2-a r (2p-l)] Ua < v < v . c 
I v~ p l 2 i-[5p/2-a r (2p+l)+a i ] < v 


(A3) 


(A4) 


Case 2: u m v ic < v a 

In Case 2, the synchrotron self-absorption frequency is above the characteristic inverse Compton frequency, v ic <C v a - 
Since inverse Compton effect dominates around the epoch of peak bolometric brightness this situation is likely to occur 
around the peak optical brightness of the SN. The synchrotron spectral peak would still be at v a and the scalings are 
given by 


!*««/.« v ic oc t 4a -- aL - 3 , 

The temporal and frequency dependence of flux /„ at observed frequency v is then generalized by 


( 1/5/2 i 4a r -a i ~l/2 
f v OC < V 5 / 2 i Aa r-OlL- 


Vm < V < 


1 /2 


(AS) 


(A6) 


INVERSE COMPTON SCATTERING OF PHOTOSPHERIC EMISSION 


The same relativistic electrons which emit synchrotron emission are also responsible for the up-scattering of pho- 
tosplreric emission dominated by IR/optical/UV photons. If the scattering electrons have the power law energy 
distribution (A r (y e ) oc 7 f T p ), the up-scattered IC radiation carries the similar form The situation under 

consider ation, IC scattering of Plan ck distributed photons by relativistic power-law electrons, has be en considered in 
detail by|Rybicki & Lightman||1979| Following equation 7.31 and 6.36 of|Rybicki & Lightman| (|1979|) we get a simple 
relation 

_ p—1 P +1 


F 


is,syn 

Pv.X 


VSyn 

vx 


Ub 

P+5 

U p £ 


x £(p) 


(Bl) 


where we have collected the natural constants h, c, q , m e and dependence on parameter p inside the function t;(p). We 
have plotted this function in Figure [8] The frequency ux is the frequency of the up-scattered photons. We use the 
subscript X to distinguish it from the inverse Compton cooling break Vi C discussed in the previous section. 
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Fig. 8.— The function £(p). 


For p = 3.0, Equation |B 1 1 reduces to the well known result 

Vsyn X Pv^syn 


V X X P„ x 


Ub_ 

Uph 


(B2) 


where factor £(p) is often n egle cted. 

The essence of Equation |B1| is that the knowledge of electron distribution and seed photon field, from radio and 
optical observations r espe ctively, is sufficient to describe the inverse Compton emission completely. For example, it 
follows from Equation B1 that for p = 3.0 the IC flux is 


Pv,X — q 24 Fv,syn 


Vx 


-1 


Uph 
Ub 


(B3) 


where Ub = B 2 /8it and U p h = L\, 0 i(t)/Airr 2 c. The ratio of energy densities Uph/Us oc L{t) in the simple case of freely 
expanding shockwave ( a r = 1), since B oc 1/r in the wind medium. Optically thin synchrotron flux in a given band 


Pu.syn t 


1 . This gives the scaling of inverse Compton flux as 

F„ t x oc L{t) x t _1 


similar to that given by Bjornsson & Fransson (2004); |Chevalier, Fransson, & Nymark (2006). From Equation B1 
can now show that for other values of p, 


(B4) 

we 


^ \ P + 5 _ P + 1 

F„ tX oc L(t) 8 xt 4 


(B5) 













